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Abstract 

We report the first observation of the radiative hyperonic B decay B + —> pAy, using a 140 fb' 1 
data sample recorded on the Y(4S) resonance with the Belle detector at the KEKB asymmetric 
energy e + e _ collider. The measured branching fraction is B{B + —> pAy) = (2.16lg'^±0.20) x 10 -6 . 
We examine its Af^- distribution and observe a peak near threshold. This feature is expected by 
the short-distance b —* sj transition. A search for B + —* 7 yields no significant signal and we 
set a 90% confidence-level upper limit on the branching fraction of B(B + —> 7 ) < 4.6 x 10~ 6 . 

PACS numbers: 13.40.Hq, 14.40.Nd, 14.20.Dh, 14.20.Jn 
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The b —> sy penguin diagram is responsible for the large rates of the observed radiative 


B —> K *7 jl] decays. It is also a good probe of new physics beyond the Standard Model 
Recently, the Belle collaboration reported a very stringent limit of 0( 10 -6 ) on the branching 
fraction of two-body B + 
body decay B° —*• pAtt~ 




I )A decays [3| but found an unexpectedly large rate for the three- 
, which proceeds, presumably, via the b —> s penguin process. 
One interesting feature of the B° —> pAn~ decay is that the observed proton-A mass M p y 
spectrum peaks near threshold. Naively, a suppression of 0(q>em) is expected for the B + — > 
pAy decay relative to B + —> pA if the former process is bremsstrahlung-like. In contrast, a 
short-distance b —> sy contribution can lead naturally to a non-bremsstrahlung-like energetic 
photon spectrum and an enhancement of M p y at low mass; the former distribution can be 
compared to the recently measured b —> sy inclusive photon energy spectrum 




These 


features motivate our study of B + pAy. Some theoretical predictions 16] for the branching 
fraction of B + —> pAy are at the 10 ~ 6 level, which is in the sensitivity range of the B-factories. 
We use a data sample of 152 x 10 6 BB pairs, corresponding to an integrated luminosity of 


140 fb 1 , collected by the Belle detector at the KERB 0 asymmetric energy e + e collider. 
The Belle detector is a large-solid-angle magnetic spectrometer that consists of a three-layer 
silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel 
threshold Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation 
counters (TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals 
located inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect K [) r mesons and to identify 

n 

muons. The detector is described in detail elsewhere |2[. 

To identify the charged tracks, the proton ( L p ), pion ( L n ) and kaon ( L K ) likelihoods are 
determined from the information obtained by the hadron identification system (CDC, ACC 
and TOF). Prompt proton candidates must satisfy the requirements of L p /(L p + L k ) > 0.6 
and L p /(L p + L n ) > 0.6, and not be associated with the decay of a A baryon. The proton 
selection efficiency is about 84% ( 88 % for p and 80% for p) for particles with momenta at 2 
GeV/c, and the fake rate is about 10% for kaons and 3% for pions. 

The prompt proton candidates are also required to satisfy track quality criteria based on 
track impact parameters relative to the interaction point (IP). The deviations from the IP 
position are required to be within 0.3 cm in the transverse ( x-y ) plane, and within ±3 cm in 
the z direction, where the z axis is opposite the direction of the positron beam. Candidate 
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A baryons are reconstructed from two oppositely charged tracks, one treated as a proton 
and the other as a pion, and must have a mass within 5a of the nominal A mass, as well as 
a displaced vertex and flight direction consistent with a A originating from the interaction 
point. To reduce background, a L p /(L p ± L n ) > 0.6 requirement is applied to the proton-like 
track. Photon candidates are selected from the neutral clusters within the barrel ECL (with 
polar angle between 33° and 128°) having energy greater than 500 MeV. We discard any 
photon candidate if the mass, in combination with any other photon above 30 (200) MeV, 
is within ±18 (±32) MeV/c 2 of the nominal mass of the n 0 (rj) meson. The above selection 
criteria are optimized using Monte Carlo (MC) simulated event samples. 


0 , 


Candidate B mesons are formed by combining a proton with a A and a photon 
each defined using the above criteria, and requiring the beam-energy constrained mass, 


A/bc = y-Ebeam — p Bl an d the energy difference, A E = E B — E heam , to lie in the ranges 5.2 
GeV/c 2 < M bc < 5.29 GeV/c 2 and —0.2 GeV < A E < 0.5 GeV. Here, p B and E B refer to 
the momentum and energy, respectively, of the reconstructed B meson, and -Ebeam refers to 
the beam energy, all in the T(4S) rest frame. Because of the A E > — 0.2 GeV requirement, 
background from B feed-down is negligible except that from B + —> pE °7 decay where E 
subsequently decays to Ay almost 100% of the time. The pE °7 events can form a nearby 
peak (shifted about -100 MeV in A E) with respect to the signal peak in the M bc — A E 
region. 

The dominant background for B + —> pAy decay is from continuum e + e _ —> qq processes, 
where q = u,d, s, c. The continuum background is evaluated with an MC sample of 120 
million continuum events. In the T(4S) rest frame, continuum events are jet-like while BB 


events are spherical. We follow the scheme defined in Ref. 
variables to form a Fisher discriminant 


• □ 


and combine seven shape 


jll | in order to maximize the distinction between 
continuum processes and signal. The variables used have almost no correlation with Mb c 
and A E. Probability density functions (PDFs) for the Fisher discriminant and the cosine 
of the angle between the B flight direction and the beam direction in the T(4S) frame are 
combined to form the signal (background) likelihood C s (£&). We require the likelihood ratio 
1Z = C s /(C s +Cb) to be greater than 0.75; this suppresses about 86 % of the background while 
retaining 78% of the signal. The optimal selection requirement is determined by maximizing 
N s /\/N s ± Nb, where N s and Nj, denote the expected number of signal and background 
events; here a signal branching fraction of 4 x 10 ~ 6 is assumed. 
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FIG. 1: The distributions of A E (for M bc > 5.27GeV/c 2 ) and M bc (for —0.135 GeV < A E < 
0.074 GeV) for B° —> pA y candidates having M p -^ < 2.4GeV/c 2 . The solid, light dashed and dark 
dashed lines represent the combined fit result, fitted background and fitted signal, respectively. 
The dotted lines represent projections of 10 assumed 7 events for comparison. 

We perform an unbinned extended maximum likelihood fit to the events with —0.2 GeV< 
A E < 0.5 GeV and M bc > 5.2 GeV/c 2 in order to determine the signal yield, £ feed-down, 
and qq background. The extended likelihood function is defined as 

~(N A +N s +N q<j ) N 

C =-- n kAP A (M bc< , A Ei) 

iV ' i =1 

+ Vy;Ps(Tfbci) AF/j) + N qq P q q(Mb Ci , AEi) , 

where N is the total number of events in the fit; P\, Pz, and P qq are the PDFs for pAy, 
pjf 7 , and continuum background, respectively; N\, and N qq are the corresponding 
number of candidates. 

The pA'y and pjf 7 PDFs are two-dimensional functions approximated by smooth his¬ 
tograms from MC simulation. We use the parametrization first suggested by the ARGUS 
collaboration Q, /(M bc ) oc M hc yjl - (A/ bc /£ , beam) 2 exp[-^(l - (M hc /E hea , m ) 2 )}, to model 
the background M bc distribution, and a quadratic polynomial for the background A E shape. 
We perform a two-dimensional unbinned fit to the A E vs M bc distribution, with the sig¬ 
nal and background normalizations as well as the continuum background shape parameters 
allowed to float. 

The A E distribution (with M bc > 5.27 GeV/c 2 ) and the M bc distribution (with —0.135 
GeV< A E < 0.074 GeV) for the region M p -^ < 2.4 GeV/c 2 are shown in Fig. [T] along with 
the projections of the fit. The two-dimensional unbinned fit gives a B + —> pAy signal yield 
of 34.llg;g with a statistical significance of 8.6 standard deviations and a B + —> p£° 7 yield 
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of 0.0 ± 4.7. The significance is defined as ^/—21n(L 0 /T max ), where L 0 and L max are the 
likelihood values returned by the fit with signal yield fixed at zero and its best fit value, 
respectively. 

We measure the differential branching fraction of pAy by fitting the yield in bins of M pA , 
as shown in Fig. 0 and correcting for the corresponding detection efficiency as determined 
from a large MC sample of events distributed uniformly in phase space. The results of the 
fits along with the efficiencies and the partial branching fractions are given in Table [I] In 
these fits, the signal yields are constrained to be non-negative. The yield is consistent with 
null signal for higher M pA bins if the non-negative constraint is removed. The observed 
mass distribution in Fig. El peaks at low pA mass, a feature seen also in B° —> pAtc~ and 
B + — > ppK + decays [f il3l ]. 



2 2.25 2.5 2.75 3 3.25 3.5 3.75 4 

M pA (GeV/c 2 ) 


FIG. 2: The differential yield for B° —> pAy as a function of M -r. 


We also study the angular distribution of the proton in the baryon pair system. The angle 
Ox is measured between the proton direction and the y direction in the baryon pair rest 
frame. Figure El shows the efficiency corrected B yield in bins of cos Ox- This distribution 


supports the b 

N cos q —N ( 

direction of the photon. We define the angular asymmetry as A = x+ 


s y fragmentation picture where the A tends to emerge opposite the 

cos 9 \ 


N cos e x+ +N cl 


; Gy, 


-, where 


N C osGx+ an d N cosdx _ stand for the efficiency corrected B yield with cos#x > 0 and cos^x < 
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TABLE I: The event yield, efficiency, and branching fraction (£>) for each bin. 


M pX (GeV/c 2 ) 

Signal Yield 

Efficiency (%) 

B (10“ 6 ) 

< 2.2 

22.7±1;| 

10.6 

1 41+0.40 
i ' 4i -0.36 

2.2 - 2.4 


9.8 

0-7412;! 

2.4-2.6 

o-o±l:| 

9.3 

0 .0012;! 

2.6 - 2.8 

o.olS:i 

9.9 

o.ooi2;22 

2.8 - 3.4 

o.ol2:l 

9.6 

o.ooi2;l 

3.4 - 4.0 

o.oll:l 

9.6 

o.ool2;! 

Total 

33.8l|;? 

- 

2 .i6l2:l 


0, respectively. The measured value for A is 0.3612;!. 
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FIG. 3: Efficiency corrected yield versus cos 9x in the baryon pair system. 

The systematic uncertainty in particle selection is studied using high statistics control 
samples. Proton identification is studied with aA-> pn~ sample. The tracking efficiency is 
studied with a D* sample, using both full and partial reconstruction. Based on these studies, 
we sum the correlated errors linearly and assign a 4.1% error for proton identification and 
4.9% for the tracking efficiency. 











For A reconstruction, we have an additional uncertainty of 2.5% on the efficiency for 
off-IP track reconstruction, determined from the difference of A proper time distributions 
for data and MC simulation. There is also a 1.2% error associated with the A mass selection 


and a 0.5% error for the A vertex selection [3]. Summing the errors for A reconstruction, we 
obtain a systematic error of 2 . 8 %. 

The 2.2% uncertainty for the photon detection is determined from radiative Bhabha 
events. For the 7 r° and p vetoes, we compare the fit results with and without the vetoes; 
the difference in the branching fraction is 0.5%, which is taken as the associated systematic 
error. 

Continuum suppression is studied by varying the selection criteria on 1Z in the interval 
0 - 0.9 to see if there is any systematic trend in the signal fit yield. We quote a 2.5% error 
for this. 

The systematic uncertainty from fitting is 2.2%, which is determined by assuming uncor¬ 
related Mbc and A E PDFs, and by varying the parameters of the signal and background 
PDFs by ±lc. The MC statistical uncertainty and modeling with six M p bins contributes 
a 4.4% error (obtained by changing the bin size). The error on the number of total BB 
pairs is 0.5%. The error from the sub-decay branching fraction of A —> p7i~ is 0.8% Q- 

We combine the above uncorrelated errors in quadrature. The total systematic error is 
9.2%. 

We see no evidence for the decay B + —> pE° 7 . We use the fit results to estimate the 
expected background, and compare this with the observed number of events in the pE°y 
signal region (—0.20 GeV< A E < 0.04 GeV and M bc > 5.27 GeV/c 2 ) in order to set an 
upper limit on the yield [15, Q, 17]- The estimated background for < 4.0 GeV/c 2 is 
84.0 ± 9.2, the number of observed events is 96, and the systematic uncertainty is 9.2%; 
from these, the upper limit yield is 35.5 at 90% confidence level. Assuming the B° —> pE °7 
three-body decay is uniform in phase space, the overall efficiency including the loss from 
the M p j < 4.0 GeV/c 2 requirement is 5.1%; the 90% confidence-level upper limit for the 
branching fraction is B(B° —> pE° 7 ) < 4.6 x 10” 6 . 

In summary, we have performed a search for the radiative baryonic decays B + 


pAy, 


^0 


and pE 7 with 152 million BB events. A clear signal is seen in the pAy mode, and we 


1-6 


measure a branching fraction of B(B + —> pAy) = (2.16 ^^53 (stat) ±0.20 (syst)) x 10“ 
which is consistent with the upper limit set by CLEO 3]. The yield of the B° —> pE°y 
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mode is not statistically significant, and we set the 90% confidence level upper limit of 
B{B° -> pT° 7 ) < 4.6 x 10~ 6 . 
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